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Abstract 

We present approximations to tree level multiparton scattering amplitudes which 
are appropriate when two partons are unresolved. These approximations are required 
for the analytic isolation of infrared singularities of n + 2 parton scattering processes 
contributing to the next-to-next-to-leading order corrections to n jet cross sections. In 
each case the colour ordered matrix elements factorise and yield a function containing 
the singular factors multiplying the n parton amplitudes. When the unresolved par- 
ticles are colour unconnected, the approximations are simple products of the familar 
eikonal and Altarelli-Parisi splitting functions used to describe single unresolved emis- 
sion. However, when the unresolved particles are colour connected the factorisation is 
more complicated and we introduce new and general functions to describe the triple 
collinear and soft /collinear limits in addition to the known double soft gluon limits of 
Berends and Giele. As expected the triple collinear splitting functions obey an = 1 
SUSY identity. To illustrate the use of these double unresolved approximations, we 
have examined the singular limits of the tree level matrix elements for e"^e~ 5 par- 
tons when only three partons are resolved. When integrated over the unresolved regions 
of phase space, these expressions will be of use in evaluating the 0{a^^) corrections to 
the three jet rate in electron positron annihilation. 



1 Introduction 



Electron-positron annihilation has proved a very clean and direct source of information about 
the nature of the strong interaction. In particular, three jet events observed at DESY gave 
the first clear indications of the existence of the gluon while more recently, the Casimirs 
that determine the gauge group of QCD have been measured using four jet data. On 
the theoretical side, a large amount of effort has been devoted to determining multi-jet rates 
within perturbative QCD. Leading order predictions for the production of up to five jets have 
been known for some time ^ |^, |^, 0] . While the general features of events are generally well 
described by leading order estimates, significant improvement can be obtained by including 
higher order corrections at either fixed order or via resummation of the dominant logarithms 
or a mixture of both. The next-to-leading order corrections to three-jet like quantities were 
first computed in the early 1980 's f^, ^ and systematically evaluated by Kunszt and Nason 
using a general purpose Monte Carlo program. In many cases the radiative corrections are 
significant and resummations of infrared logarithms must be employed to make sensible 
comparisons with experimental data. Such calculations have been used to extract a precise 
value of the strong coupling constant from three jet events and other hadronic event shapes 
T0| with a global average of. 



as{Mz) = 0.121 ±0.005. 
A recent re-evaluation of lower energy data yields ||TT[| , 



as{Mz) = 0.122+°-°°« 



'-0.006; 



while preliminary analyses using hadronic data above the Z-resonance from the OPAL Col- 
laboration |[T^ shows that the strong coupling evolves in the manner predicted by QCD. 

Recently there has been much progress towards a more complete theoretical description 
of four jet events produced in electron-positron annihilation. The one-loop matrix elements 



relevant for four jet production have been evaluated by two groups |]T3|, |14|, |T5|, |16| using 



contrasting methods. Bern et al calculate helicity amplitudes for the e~^e (7*,^) 



4 parton processes while in Refs. [13, 15|, the spin summed interferences between one-loop 



and tree-level 7* — 4 parton processes were evaluated. These results have been compared 



T6| and found to agree for specific phase space points. Using these one-loop matrix elements 



and the earlier tree level five parton matrix elements of [0, two general purpose Monte Carlo 



programs have been constructed ||T^, |T8[ and used to compute the four-jet rate and other 
four-jet like hadronic observables at next-to-leading order (i.e. 0{al)). The corrections 
are again large |T^, [T^ indicating the need for resummations of infrared logarithms, but 



nevertheless, improved determinations of the Casimirs are possible [19| together with more 



stringent bounds on the existence of a light supersymmetric partner of the gluon pU|. 

One of the next steps is to evaluate the next-to-next-to-leading order corrections to the 
three-jet rate. Together with LEP/SLC data and hadronic events from the next linear 



1 



collider this would play a role in reducing the error on as{Mz) from event shape analyses 
to the 2-3% level. However, to accomplish this, several ingredients are necessary. First, 
the two-loop 7* — > 3 parton amplitudes must be determined including the evaluation of 
the two loop box graph with massless internal and external legs. Second, the real infrared 
singularities which occur when two particles are unresolved in the tree level process and 
when a single particle is unresolved in the one-loop process must be systematically isolated. 
In the case of two real unresolved particles, there are a variety of different configurations; 



(a) two soft particles, 

(b) two pairs of coUinear particles, 

(c) three collinear particles, 

(d) one soft and two collinear. 



while in the one-loop process one particle may be soft or two may be collinear as before. 
For next-to-leading order calculations, many methods have been developed to first isolate 
the infrared singularities associated with one unresolved particle analytically and then nu- 
merically combine the remaining finite real and virtual contributions 0, ^ ^, |2^, 23]. All 



of these rely on suitable approximations to the real matrix elements in the soft and collinear 
limits. However, the infrared singular limits of the real and one-loop amplitudes relevant 
at next-to-next-to- leading order are rather less well known. In ref. [^, Bern et al., have 



developed appropriate splitting functions for one-loop processes where two external particles 
are collinear while Berends and Giele have examined the multiple soft gluon behaviour of 
QCD matrix elements in ref. ||2^ . 



In this paper, and as a first step towards the next-to-next-to-leading order 3 jet rate, 
we examine the double unresolved limits of multiparton scattering amplitudes, and find 
suitable approximations for the e'^e" —>■ 5 parton matrix elements when two particles are 
unresolved. We present results for all of the configurations (a)-(d) above. In section ^ we 
briefly review the structure of tree level scattering amplitudes when one particle is unresolved. 
In section |] we write down general expressions for the structure of the e~^e~ — >■ 5 parton 
matrix elements. Sections ^ and ^ are organised according to whether the two unresolved 
particles are colour connected or not. The precise meaning of what colour connected means 
will be given in sect. ^. In the unconnected case, the singular limits are merely obtained 
by multiplying single unresolved factors. However, when the particles are colour connected, 
the structure is more involved (sect. |^) and we give explicit formulae detailing the double 
unresolved singular factors. In each case we write down expressions for the double unresolved 
limits of the five parton process. Our results are summarised in section M. 
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2 Single unresolved limits 



The soft gluon and coUinear parton limits of multiparton scattering amplitudes are well 
known. For example, if we consider the matrix element for the tree level scattering of n 
gluons labelled 1 . . . n, then, decomposing the amplitude according to the various colour 
structures [EBI, we find, 



>i(l,2,...,n) = 2i^"-2 TriT'^'T"^ ...T^")A{l,2,...,n). (2.1) 

P(l,2,...,n-1) 

Here, represents the colour of the ith gluon and the summation runs over the (n — 1)\ 
non-cyclic permutations of the gluons. At leading order in the number of colours, N, the 
squared matrix elements summed over all colours have the simple form, 

(2 ]\r\ 
V ^^'-1) E 1^(1,2,. ..,n)|l (2.2) 
^ / P(l,2,...,n-1) 

Unlike QED , it is the colour ordered sub-amplitudes A rather than the squared matrix 
elements that have nice factorisation properties in the soft and coUinear limits (see for 
example pSf). For example, in the limit where the nth gluon becomes soft, we have the 
QED-like factorisation into an eikonal factor multiplied by the colour ordered amplitude 
with gluon n removed, but the ordering of the hard gluons preserved, 

1^(1,2, ... ,n - l,n)p 5'(„_i)„i(si(„_i),si„,s„(„_i)) \A{1,2, . . . ,n - 1)|^. (2.3) 

with the eikonal factor given by, 

Sabc{Sac,Sab,Sbc) = —■ (2-4) 

^ab^bc 

Similarly, in the limit where two gluons become coUinear, the sub-amplitudes factorise. For 
example, if gluons a and b become coUinear and form gluon c, then adjacent gluons give a 
singular contribution, 

1^(1, . . . , a, 6, . . . , n) |2 ^ Pgg^giz, Sab) \A{1, . . . , c, . . . , n)\^ , (2.5) 

while separated gluons do not, 

1^(1,..., a,..., 6,..., n)|2 ^finite. (2.6) 

In other words, there is no singular contribution as Sab 0, and, when integrated over the 
small region of phase space where the coUinear approximation is valid, this colour ordered 
amplitudes generates a negligible contribution to the cross section. In eq. (2.5), z is the 
fraction of the momentum carried by one of the gluons and, after integrating over the az- 
imuthal angle of the plane containing the coUinear particles with respect to the rest of the 
hard process, the coUinear sphtting function Pgg^g is given by, 

P99^gi^,s) = -^P„{z) (2.7) 
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where the usual Altarelh-Parisi sphtting kernel with the colour factor removed is given 
by, 

In either soft or collinear limits, the process appears to involve only n — 1 hard gluons and 
the infrared singularities must precisely cancel 0, |31|] with those from the [n — l)-gluon 
one-loop amplitudes. Together, the real and virtual graphs form the next-to-leading order 
perturbative contribution for infrared safe quantities associated with the {n — 1) gluon scat- 
tering process. Many methods have been developed to first isolate the infrared singularities 
analytically and then numerically combine the remaining finite real and virtual contributions 

i, in, I, nil. 



3 Colour-ordered matrix elements squared 

There are two five parton processes, 

e+e" qqggg, (3.1) 

and, 

e+e~ qqQQg, (3.2) 
for which the lowest order matrix elements for process ( |3.1|) can be written, 



M{Qi; 1, 2, 3; Q,) = SUQ^, 1, 2, 3; Q,)V^, (3.3) 



while for process ( p.2|) we have, 

M{Qi,Q,; Qs,Q,; 1) = f.'iQuQ^; Q3,Q,; 1)^^ (3.4) 

In these expressions, represents the lepton current, while and are hadronic currents 
containing quarks and gluons. The gluon colour is denoted by the adjoint representation 
label oi, . . . , a3 while that of the quark is ci, . . . , C4. 

The two quark-three gluon current 5^ may be decomposed according to the colour struc- 
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SliQu 1, 2, 3; Q2) = ieg"" ^ (T"^T''^T'^^),,,,5^(Qi; 1, 2, 3; Q^), (3.5) 

P(l,2,3) 

where S^{Qi; 1, 2, 3; Q2) represents the colourless subamplitude where the gluons are emitted 
in an ordered way from the quark line. The colour matrices are normalised such that, 

Tr (T"r^) =^S''^- 



By summing over all permutations of gluon emission, all Feynman diagrams and colour 
structures are accounted for. 

The squared matrix elements for eq. are simply, 



X 



N 



E 

P(l,2,3) 
+ 



S^{Qi;l,2,3;Q,)V>^ 
+ 1 



S,{Q,-l,2,~3-Q,)V'^ 



S^iQ,; 1,2,3; Q,)V^ 



(3.6) 



In the last two terms, the tilde indicates that that gluon should be inserted in all positions 
in the amplitude. In other words. 



S^iQu 1, 2, 3; Q2) = MQi; 1, 2, 3; Q^) + 5^(Qi; 1, 3, 2; Q^) + S^iQi; 3, 1, 2; Q, 



(3.7) 



In this case, gluon 3 is effectively photon-like and the contribution from the triple and quartic 
gluon vertices drops out. For reference, the squared matrix elements for quark- ant iquark 
production together with a single gluon can be written. 



2 ^^2 92L\ (N'~l 



N 



(3.8) 



The four quark-one gluon current may be decomposed according to its colour structure 
as follows, 

3 



X 



t:iJ,,,MQi, Q,; Qs, Q,; 1) + (Qi ^ Qs, Q, ^ Q,) - {Qi ^ Q3) - (Q^ ^ Q,) 



(3.9) 

where the exchanges are understood to apply to the colour labels as well. The colour ordered 
sub-current can be written, 



^(Qi! Q2; Qs, Q4; 1) — '^QiQ2'^Q3Q4^j (Qi) Q2; Q4; 1) 

+ ^'^Q3(32^giQ4^^(Qi) Q2; Q3, <34; 1)) 



(3.10) 



where, 



r/(gi,Q2;Q3,Q4;i) = ^^'^HQi; i;Q4lQ3;Q2) +^?^^(Q3;Q2lQi; i;Q4) 
T/(gi,g2;Q3,g4;i) = S?^^^(gi;i;g4lQ3;g2) + S?^'(Q3;g2|gi;i;Q4)- (3.11) 



Here, 5q^q^ = 1 if quarks i and j have the same flavour. The functions A^^^'^ and B^^^'^ 
describe Feynman diagrams where the gauge boson couples to the Q1Q2 pair. However, in 
A^^'^'^, the colour flows along the gluon connecting the two quark pairs, so that Qi and 
(^4 are the endpoints of a colour antenna (and similarly Q3 and Q2) while in B'j^^^'^, no 
colour is transmitted between the quark pairs and now Qi and Q2 form the endpoints of a 
colour antenna (and similarly Q3 and Q4). In each case, the gluon may be emitted from any 
position on the colour line. 

Squaring the four quark-one gluon amplitude and summing over colours yields. 



'g^NY (N^ 



X 



7V2 



%{Qi. Q2\ Qs, Q4; + MQ,, Q,; Q,, Q2; 1)V^ 



(r,{Q,, Q2; Q3, Q,; l)V>' + r,{Qs, Q,- Qi, Q2; 1)^^) [%{Qi. Q4; Q3, Q2. 



N 

+((5i ^ (53,<32 ^ Qa)-, 

or equivalently. 



(3.12) 



X 



+ 



T^^(Qi,Q2;Q3,Q4-A)V'' 
1 



rfiQi, Q,; Qs, Q2; 1)V'' ' - T,{Q,, Q^; Qs, Q,; 1)V^^ 



Q1Q3 



N 

t2 



K (T/(gi, Qs, Q4; ^)v^) {t,^{Qu Q2, Qs, Q4; m")^ 



+{Qi ^ Qs, Q2 ^ Qa) + ^QiQsiQi ^ Qs) + ^Q2Qa{Q2 ^ Qa) 



(3.13) 



Here we have introduced the shorthand notation. 



'r^(Qi,Q2;Q3,Q4;i) = t/(Qi,Q2;Q3,Q4;i)+'?'/(Q3,Q4;Qi,Q2;i) 

= rf{Q,,Q,;Qs,Q2;l)+rf{Qs,Q2;Qi,Q,;l). (3.14) 

Note that in the case of identical quarks, there is an extra symmetry factor of 1/4 multiplying 
the matrix elements. 
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4 Colour unconnected double unresolved 



In the cases where the two unresolved particles are not colour connected, the factorisation of 
the amplitudes involves the well-known functions describing single soft and collinear emission. 
We first describe what is meant by colour connected and colour unconnected. 



4.1 Colour connection 



In the previous section we have seen how tree level matrix elements can be decomposed into 
colour ordered subamplitudes which have nice factorisation properties in the infrared limits. 
Therefore it is useful to view the matrix elements in terms of the colour structure associated 



with the subamplitudes. For example, in eq. ( |3.5| ), the two quark-three gluon subamplitude 



I5 2, 3; Q2) is associated with the colour structure {T"-^T"''^T"-^)c^c2- This is a colour 
antenna that ends on the quark/ ant iquark colour charges Ci and C2 with ordered emission of 
gluons with colour ai, . . . , 03. Within this colour antenna, gluon 1 is colour connected to the 
quark Qi and gluon 2, but not to the antiquark Q2 or to gluon 3. In cases involving more 
than one quark-antiquark pair there can be many colour antennae. For example, the four 
quark amplitude A^^'^'^{Qi; 1; Q4,\Q3', Q2) of section 3 describes a process with two separate 
colour antennae. In general, the particles in one antenna are not colour connected to the 
particles in one of the other antennae. However, there is one case where particles in adjoining 
antenna can usefully be thought of as colour connected. This is when there is an antiquark 
at the end of one antenna and a like flavour quark at the beginning of another, 

^(...,Q|Q,...)- 

When this quark-antiquark pair are collinear, they combine to form a gluon G, which then 
connects, or pinches together, the two separate colour antennae, so that, 

|^(. . . , QIQ, . . ^ P,5^g(^, Sqq)\A{. ..,a,G,b,.. .)\\ 

A useful definition of colour "connected" therefore includes these antennae pinching config- 
urations along with the more straightforward colour connection within a single antenna. All 
other cases are colour "unconnected" . 



4.2 Two collinear pairs 

Two pairs of particles may become collinear separately, but with the particles in one or 
both of the pairs themselves not colour "connected" . In these cases, there are no singular 
contributions containing both of the vanishing invariants. For instance, if partons {a, d} and 
{b, c} are collinear then, 

|^(. . . , a, . . . , 6, . . . , c, . . . , . . .)|^ — i> less singular. (4.1) 
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By this we mean there is no contribution proportional to 1/sadSbc and once again, when 
integrated over the small region of phase space relevant for this approximation yields a 
negligible contribution. 

The situation where two pairs of colour "connected" particles are coUinear is also rather 
trivial. If partons a and h form P, while c and d cluster to form Q, so that P and Q are 
themselves colour unconnected, then. 



|^(. . . ,a,6, . . . ,c,rf, . . .) 



Pab^p{Zl, Sab) Pcd^Qiz2, S,d) |^(. . . , P, . . . , Q, . . .) (4.2) 



Here, zi and Z2 are the momentum fractions carried by a and c respectively. A similar result 
holds if either of the pairs involves particles in separate antennae, but which are able to 
undergo antenna pinching. The coUinear splitting functions are related to the (colourless) 
Altarelli-Parisi splitting kernels by eq. ( |2.7|) which, in the conventional dimensional regular- 
isation scheme |^ with all particles treated in ci = 4 — 2e dimensions, are given by ||29|| , 



qg^q 



qq~*g 



'l + z"" -e{l-z) 
'z^ + jl-zf-e 



(4.3) 
(4.4) 



with Pgg_g glvcu lu eQ. (^) and Pgq^q{z) - . qg-,g 
of the coUinear particle plane is understood 



Pqg-,g{l — z). As before, azimuthal averaging 



4.2.1 Double coUinear limit of e+e — > 5 partons 

In this limit where two pairs of partons are simultaneously coUinear, the five parton matrix 
elements factorise into the three parton matrix elements multiplied by a combination of 
products of coUinear splitting functions. Summing over all possible unconnected double 
coUinear limits, for the two-quark currents we find. 



5^(Qi;l,2,3;Q2)^'' 
5^(Qi;i,2,3;Q2)V^^ 



-Pqi1^Q-P23^G + Pq^I^qP^q^^q + Pl2^GP-iQ^^- 



Pqi1^qP2Q2^q + Pqi3^qP2Q2^q + Pqi1^qPzq^->q 
+Pqi3^qPi2^g + Pi2^gPzq^^q) '5^^'' 



51 Pqi 

P(l,2,3) 



i^qP2Q2^q 



(4.5) 



whilst the only contributing pieces for the four-quark process are. 



rf{Q^,Q,;Q,,Q,;l)V^ 



^f(Qi,Q4;Q3,Q2;i)^^ 



PQi1^qPq^Q3-*G + PlQi^QPQ2Qi^G 
{PqiI^qPq^Qs^G + PlQ2^Q^Q4Q3^G 
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P — —P— 



+Pq31^qPq2Q^^g + Piq^^qPq^q.^g) l^f^l ■ (4-6) 

For brevity we have dropped the arguments of the sphtting functions. Exphcitly, for np 
flavours of quark, we find, 



3! 



+ 



1 



Ar2 



12-*G 



3Q2^Q 



2nF 

'If' 



Pt^ 



^ (Pq,i-Q^2Q,^q + Pqi3^qP2Q,^q) ) • (4-7) 



2iV2 



Note also that the identical gluon factor 1/3! is eliminated since each term in the sum over 
permutations produces an identical contribution. It is interesting to review the origin of 
the factor of np. There are np{nF — l)/2 contributions from two unlike pairs of quarks, 
each of which generates two sets of singular limits - that indicated plus the symmetric term 
{Qi ^ QZ1Q2 ^ Qa)- addition there are np like-quark pair contributions, which after 
the symmetries have been applied yield four singular limits. However, the identical quark 
contribution is multiplied by the identical particle factor 1/4 so that the net result is. 



nplnp 



1) . . 1 

— x2 + nFx4x- 

4 



One factor of np is absorbed into the three parton matrix elements 
appears as an explicit factor. 



, while the other 



4.3 Triple coUinear factorisation 

If three collinear particles are colour "unconnected" then there is no singularity. So if a, b 
and c all become collinear, 

|^(. . . , a, . . . , 6, . . . , c, . . .) 1^ ^ finite, (4.8) 

and there is no singular contribution involving the invariants Sab, Sbc or Sabc- As before, 
because the region of phase space where the triple collinear limit is valid is extremely small, 
this gives a negligible contribution to the cross section. When two of the three collinear 
particles are colour "connected" we find a singular result, 

\A{...,a,...,b,c,...)\' ^1/sbc. (4.9) 

However, when integrated over the triple collinear region of phase space that requires Sab, Sbc 
or Sabc all to be small, we again obtain a negligible contribution that is proportional to the 
small parameter defining the extent of the triple collinear phase space. We therefore ignore 
contributions of this type. 
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4.4 Soft /collinear factorisation 

Two particles may be unresolved if one of them is a soft gluon and another pair are collinear. 
When the soft gluon g is not colour connected to either of the colour "connected" collinear 
particles c and d, factorisation is straightforward, 

|^(. . . ,a,g,b, . . . ,c,d, . . .)\'^ Sagb{sab, Sag, Sbg)Pcd^p{z, Scd) \ A{. . . , a, 6, . . . , P, . . or 

(4.10) 



4.4.1 Soft/collinear limit of e^e 5 partons 

In the soft/collinear limit, the five parton matrix elements again factorise into a singular 
factor multiplying the squared two-quark current relevant for three parton production. 



5^(Qi;l,2,3;Q2)^^ 
5^(gi;l,2,3;g2)V^ 



'S'Qi12^3Q2^q + Pq 



(4.11) 



Note that for 



S,iQi;l,~2,~3;Q,)V^ 



the soft and collinear limits are considered to be over- 



lapping and will be dealt with in section ^.2| . 



In the four-quark current case, the soft/collinear limit has only two colour- unconnected 
contributions. The first is given by. 



rfiQuQ,;Q3,Q2-A)V^ 



(4.12) 



whilst the limit of, Tfj,(Qi,Q2',Q3,Q4^',l)V^ , again involves both unconnected and con- 
nected factors and therefore discussion of this will also be deferred until section 5^. The 
other subamplitudes vanish in the unconnected soft/collinear limit. 



Applying these limits to the full five parton matrix elements is straightfoward and, after 
removing identical particle factors where necessary, we find. 



1 
3! 



+ 



SqiuP; 



PQi1^qS23Q^^ 



1 

iV2 



Pi2~*gS, 



+ !!£C -p - 



(4.13) 
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4.5 Two soft gluons 



When two unconnected gluons are soft, the factorisation is again simple For gluons gi 
and g2 soft we find, 



|^(. ..,a,gi,b,...,c, g2,d, 



x|^(...,a,6,...,c,ci,...)r, (4.14) 



so that the singular factor is merely the product of two single soft gluon emission factors 
given by eq. (|2.4|). Note that 6 = c is allowed. 



4.5.1 Double soft limit of e+e 5 partons 



The sum over the unconnected double soft limits of the colour ordered subamplitudes can 
be easily read off. 



S^{Qr- 1,2,?,- Q^)V^ 
5^(Qi;l,2,3;Q2)V^'^ 
5^(gi;i,2,3;g2)^'^ 



5'qi12'S'23Q, 



9 ^QllQi^Ql 



fj. 



2Q2 



P(l,2,3) 



(4.15) 



There is no contribution from the four-quark matrix elements. Inserting these limits into 
the full five parton matrix elements yields. 



1 
3! 



+ 



TV 



2 fg^NV 



X 



(4.16) 



2N^ I "'^ilQa <9i2Q2 

where once again the sum over permutations is eliminated by the identical particle factor. 
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5 Colour connected double unresolved 



The factorisation that occurs when the two unresolved particles are colour "connected" is 
necessarily more involved than that in section ^ In particular, we will need to intoduce new 
functions to describe this factorisation. 



5.1 Triple coUinear factorisation 

When three colour "connected" particles cluster to form a single parent parton there are 
four basic clusterings, 

999 G, qgg Q, 
gqq G, qqq -> Q, 

and the colour ordered sub-amplitude squared for an n-parton process then factorises in the 
triple coUinear limit, 

|^(. . . , a, 6, c, . . .) I' ^ Pabc^p\A{. . . , P, . . .) I'- (5.1) 

As before, partons able to undergo antenna pinching are considered to be colour connected, 
so that there may be contributions from amplitudes such as A{. . . , a, b\c, . . .). The triple 
coUinear splitting function for partons a, b and c clustering to form the parent parton P is 
generically, 

Pabc-^p{w, X, y, Sab, Sac, Sbc, Sabc), (5.2) 

where w, x and y are the momentum fractions of the clustered partons. 

Pa = wpp, pb = xpp, pc = ypp, with w + X + y = 1. (5.3) 

In addition to depending on the momentum fractions carried by the clustering partons, 
the splitting function also depends on the invariant masses of parton-parton pairs and the 
invariant mass of the whole cluster. In this respect, they are different from the splitting 
functions derived in the jet-calculus approach B^, and implemented in the shower Monte 



Carlo NLLJET [Q, which depend only on the momentum fractions. 



The triple coUinear splitting functions Pabc^p are obtained by retaining terms in the full 
matrix element squared that possess two of the 'small' denominators Sab, Sac, Shc and Sabc- 
As before, we consider the explicit forms of the 7* — > four and five parton squared matrix 
elements and work in conventional dimensional regularisation, with all external particles 
in d = 4 — 2e dimensions. Similar results could be derived using helicity methods or by 



examining the on-shell limits of the recursive gluonic and quark currents of ref. . 



Although the splitting functions are universal, and apply to any process involving the 
same three colour connected particles, for processes involving spin-1 particles, there are ad- 
ditional (non-universal) azimuthal correlations due to rotations of the polarisation vectors. 
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These angular correlations do not contribute to the underlying infrared singularity structure 
and vanish after all azimuthal integrations have been carried out and we therefore system- 
atically omit them. 

A further check on our results is provided by the strong-ordered limit, where the particles 
become collinear sequentially rather than at the same time. In this limit these functions 
factorise into the product of two usual collinear splitting functions plus azimuthal terms, 



agreeing with the results of [35] 



5.1.1 Three collinear gluons 

Firstly, examining the sub-amplitudes for multiple gluon scattering, we find that the colour- 
ordered function Pggg^Q is given by, 

Pabc^ciw, X, y, Sab, Sbc, Sabc) = 8 X | 

l-e){xsabc-{l-y)sbcf , 2{1 -e)sbc , 3(1 - e) 



^ab^abc ^abS^bc 



abc 



^ 1 / (l-y(l-y))^ ^ x'^ + xy + y'^ ^ xw-x'^y-2 ^ x 



SabSabc\ yw{l-w) 1-y yi^-y) i^-y), 

^ 1 L^2 2i2-w + w')ix' + w{l-w)) ^ 1 ^ 1 



2sabSbc\ y{^-y) yw {i-y){i-w)J} 

+ {sab ^ Sbc, w y) + azimuthal terms. (5.4) 

This splitting function is symmetric under the exchange of the outer gluons, and contains 
poles only in Sab and Sbc- 



5.1.2 Two gluons with a collinear quark or antiquark 

There are two distinct splitting functions representing the clustering of two gluons and a 
quark which depend on whether or not the gluons are symmetrised over. In the unsym- 
metrised case, there will be poles in Sg^g^, due to contributions from the triple gluon vertex 
which are not present in the QED-like case. For the pure QCD splitting we find, 

-^i?9iff2— ►Q ('^5^)1/5 ^qgi 5 '^952 5 '^9192 ) ^q9i92 ) ~ ^ X 

^991^9192 v V y i^-'^) J \y ^-wjj 

+ f (1 - f (1 - + --^-^y)., f - ^)'^ - - - .^.) 

^991*99192 V V y{^-w) J \ y{^-w) J J 
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'1 — yY{2 — y) + + 2xw — 2 — y\ ^ {xw — y — 2yw)\ 

S9ig2Sqgig2 V 1/(1 -W) ; y{l-w) ) 



1-e 



+ ( '^";r;. + ^ (4^ + d - + (3 - ^)) ) U-^) 

while for the QED-hke sphtting where one or other or both gluons in the colour ordered 
amplitude are symmetrised over, 

"^99192^0 (^? •^5 Hi ^qgii ^qg2i ^qgig2^ 4 X 

1 W / o / 2 2\ 2 \ 

+ + w — e[x +xy + y)—exyj 

'^^qgi^qg2 -^V 

+ (""^(1 — X + e^xy) + (1 — y)'^ — e(l — y){x^ + xy + y"^) + e^xy^ 

'^^ '1 - e)^ - e ) I + ^ ^ y). (5.6) 



^qgig2 \ ■^992 



The function P^g^g^^q can be interpreted as the relevant triple coUinear splitting function 
with one or both of the gluons replaced by photons. As such, this result echoes that found 
in for Pqg^^q. Using charge conjugation, we see that the functions representing clustering 



of two gluons with an antiquark are simply, 

^gig2q^Q^^ ^ ^' "^91925 -^9295 ^giq^ -^91929) ~ Pqgigi^QiVi ^5 ^? -^5925 ■^$91' '^91925 -^91925)5 

Pg^g-iq^qi'^^-^^y^ ^g2q^ ^giq^ ^gig2q) ~ Pqgig2^Qiyi-^^'^i^qg2^^qg\^^gig2q)- (S-"^) 



5.1.3 A quark-antiquark pair with a collinear gluon 

Similarly the clustering of a gluon with a quark-antiquark pair into a parent gluon again has 
two distinct functions. For example, there is a singular contribution from the four quark 
matrix elements when (^4, and the gluon cluster, 

r^{Qi.Q2;Qz.Q,-A)V^^ = ^?'^^(Qi;i;Q4lQ3;Q2)V^^ + ^?'^'(Q3;Q2lQi;i;Q4)l^' 



Here, the singular contribution comes entirely from the first A term where combining Q^, 
Q3 pinches the two colour lines together. Explicit evaluation yields. 



■^gqq^ci.^ 1 X, y, Sgg, Sgg, Sggg) 4 X 



^gqq \ ^qq 



+ (l-e)^ + (3- 



2(xs 



gii 



'1 — w)s 



gqj 



^gi 
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+ 
+ 
+ 



1 



^ ■-'^ -y-2w-e- ^ ^^^^ ' 



w{l — w) (1 — e)w{l — w) ) 

1 (x{{l-wf-w^) 2x'^ {I - yw - {I - y){l - w))" 



SgqSqq 



w(l — w) 



[l-e)w{l-w) 



(1 + + 4a;w) 2x {w{x — y) — y{l + w)) 



+ 



SqqSgqq 



w{l — w) (1 — e)w{l — w) J \ 

Again applying charge conjugation yields the further relation, 

Pqqg^ciw 1 X, y, Sqq, Sqg, Sqqg) = Pgqq^Q(jJ , X, W , Sgg, Sgq, Sqqg^, 



azimuthal terms. (5.8) 



(5.9) 



describing instances where the gluon is colour connected to the quark rather than the anti- 
quark. 

There is a further contribution when the quark- antiquark and gluon combine to form a 
photon-like colour singlet. This occurs when, 



V(^3,Q2;Qi,Q4;i)^^ 



In this case the singular contribution is produced by the second B term and is due to the 
entire Qs; 1; antenna pinching to form a gluon which is then inserted in a symmetrised 
way (i.e. with a tilde) into the remaining colour antenna. This QED-like splitting function 
is given by. 



1 



:i-e) 
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+ 1 + 



SgqSqg 



^qg^qgq 



\ + 2x + e y-^ I -|- {sqg Sgq, w ^ y) + azimuthal terms.(5.10) 



5.1.4 A queirk-antiquEirk pair with a collineEir queirk or antiquark 

Lastly, we consider the clustering of a quark-antiquark pair (QQ) and a quark (q) to form 
a parent quark q' with the same flavour as q. The splitting function depends upon whether 
or not the quarks are identical, 

p pnon— idcnt. qQ pidcnt. (^^^\ 

UQQ^q' - -^qQQ^q' 'JfqQQ^q" ^^'^^^ 
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where SqQ = 1 for identical quarks. If quarks Qi, Q3 and (^4 are clustered to form Q, then 
we find a non-identical quark contribution, 



^non— ident. 
Q1Q4Q3— >Q 



with, 



- 

9QQ 



2s 



qQ 



s^-s^ -(1 
QQ qQQ^ 



+ 



■^QQ^qQQ 



+ + + 



-efl - w) 



(5.12) 



The singular contribution is now generated by the square of the second A term; pinching 
and together connects the two colour antennae together and combining with Qi ensures 
that the vector boson couples to a flavour singlet pair. Precisely the same function 

describes the triple coUinear hmit of the functions. We find that in the same hmit, 

T^iQ.^Q^^Q^^Q^AW' = Bj^«^(Qi;i;g2lQ3;Q4)v^'' + ^?^^(Q3;Q4lQi;i;Q2)v^'^ 



^non— ident. 
QiQ^Qs^Q 



where this time the first term alone contributes. Here Q3 and combine to form a photon 
which chistcrs with Qi. As before, charge conjugation generates the associated function for 
an antiquark combining with a quark- antiquark pair. 



-P^Q^_l^^''*'('f^> ^' ^' *9<3' *QQ' ^gQg) ~ ^q^Q-^q'^^'('^^ ^' ^' ^qQ^ ^QQ^ ^qQQ)- 



(5.13) 



When the flavours of the clustering quarks are the same, there is an additional contri- 
bution coming from the interference terms of the four-quark matrix elements. For instance, 
when Qi, (^4 and combine, 

- ^?'^^(Q3;i;g2lQi;Q4)^'^(^?'^^(Q3;i;g2lQi;Q4)^'^)^ 



where. 



jjident. 

2 QiQ^Qs^Q 



[jident. 
qQQ^q' 



,,{w, X, y, s^Q, Sqq, s^qq) — 4 X 
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e 



+ (^5gQ ^ ^ w^)- (5-14) 

Here, there are poles in the matrix elements when clusters with both Q3 and Q\ and the 
triple coUinear function is symmetric under q ^ Q. 



5.1.5 The N = 1 SUSY identity 

These triple-collinear splitting functions, like the ordinary Altarelli-Parisi splitting kernels, 
can be related by means of an = 1 supersymmetry identity. In unbroken supersymmetric 
theories, the masses of gluon and gluino are identical thereby ensuring that the self-energies 
of the two particles are equal. Considering all two particle cuts of the one-loop diagrams 
contributing to these self-energies, then in terms of the colour stripped Altarelli-Parisi kernels 



eq. |2]3 and eq. O 



Pgg^ai^) + Pm^ai^) = Paa^^i^) + Pa^-^M- (5.15) 

Here the quark plays the role of the gluino. Note that in conventional dimensional regular- 
isation the number of degrees of freedom for the gluon and gluino are no longer equal and 
the supersymmetry is broken^. Therefore this result is not true away from four dimensions. 
Similarly, by considering the three particle cuts of the relevant two-loop contributions to the 
self energies, and omitting the invariants in the arguments of the functions, we find. 



J2 Paaa^cia, b, c) + 2Pgqg^G{a, b, c) + Pqgq^cia, b,c) \ = 

{a,b,c) \ / 

Y: (2P,,,^Q{a, b, c) + Pqgg^Q{a, b, c) + 2P-"_;g-Ma, c) + P^'^^^Qia, b, c) 

(a,b,c) \ 



(5.16) 



P(a,6,c) 

provided we set e = 0. This non-trivial relation between the splitting functions is a further 
check on the results presented in this section. 



5.1.6 Triple collinear limit of e^e —>■ 5 partons 

Examining the two-quark currents in the triple collinear limit is now straightforward. Sum- 
ming over all triple collinear limits, we find. 



S^iQ,- 1,2,3- Q^)V' 



Pc 



P^ 



P23Q2^Q 



^We note that there are other variants of dimensional regularisation where the gluon and gluino degrees 
of freedom are equal. 



17 



S,{Q^,1,2,~3;Q,)V^ 
S,{Qi;l,~2,~3;Q,)V>' 



(-PqiI; 

9 X] (-PQii2-^Q + -^23^2-^^) 



P{1,2,3) 



(5.17) 



The four-quark current contributions are more complicated, but in each case yield factors 
multiplying the two- quark current, 



rf{Q,,Q,;Qs,Q,;l)V^ 



_j_ pnon— ident. i pnon— ident. \ 



2V3V4— 'V V3V2'' 

(p - -I- pnon— ident. _i_ pnon— ident. 1 

non— ident. i Dnon— ident. 

2^ 



p _ I pnon— lacnt. _i_ pn 

_l_ p _|_ pnon— ident. i pnon— ident. \ 



hQiQ2-*Q Q3Q2 

^ / pident. i pident. 1 pident._ i_ pident._ \ 

2 V Qi<32Q3-*Q QzQiQi^Q Q2Q3Qi-^Q QiQiQ2-^Q) 



^ / pident. 

V QiQ2< 



I pident. 

2 V" WiQ2Q3^Q Q2QzQi-^Q 



(5.18) 



Combining these limits and eliminating the identical particle factors where appropriate pro- 
vides the triple coUinear singular factor for the five parton squared matrix elements. 



1 
3! 



+ 



X 



jy2 I V '^il^^'S ^ 12Q2-*Q 



+!^(p- +p- -^p - - 

N \f^Q4Q3^G ^ ^Q^Qsl^G j^2^Q3lQ4^G 



^ N [ m 



1 



pnon— ident. _|_ pident. _i_ pnon— lacnt. piacnt. /tr 

QiQiQi^Q 27V *5iQ4Q3^Q ^ QiQiQi-^Q ~ 2iV^2Q3Q4^Q^ 



jnon— ident. 



1 



lident. 
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5.2 Soft/coUinear factorisation 

Wc now examine the configurations where one gluon is soft and two particles are coUinear. 
In this case, the sub-amplitudes factorise as. 



|^(. . . , d, . . . , a, 6, C, . . .) I' ^ Pd;ahc\A{. . . , d, . . . , P, . . .) I', 



(5.20) 
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where gluon a is soft, partons b and c are coUinear and either colour connected or able to 
undergo antenna pinching. Parton d is the adjacent colour-connected hard parton in the 
antenna containing the soft gluon. If a is symmetrised over (a, so QED-like) then d is the 
quark in that colour-line; otherwise d is simply the parton adjacent to a. In this limit the 
collinear partons form parton P and carry momentum fractions, 

Pb = xpp, pc = (1 - x)pp, (5.21) 
and we write the soft /collinear factor as, 

Scd). (5.22) 

To determine the limiting behaviour, we again consider the explicit forms of the 7* — *• 
four and five parton squared matrix elements. All terms that possess three of the 'small' 
denominators Sab, Sad, Sbc and Sabc contribute in the soft /collinear limit. Similar results 
could be derived using helicity methods or by examining the on-shell limits of the recursive 
gluonic and quark currents of ref. [^. Alternatively, these limits can be obtained directly 
from the triple collinear limits of sect. 5.1, by keeping only terms proportional to 1/w and 
subsequently replacing 1/w by {sbd + Scd)/sad, 1/ (1 — w) by 1 and y hj 1 — x. 

In fact, in this limit we find a universal soft factor multiplied by a collinear splitting 
function, 

Pd;abcix^ ^aby ^bcy ^abcy ^ady ^bdj ^cd) Sd-abcix, Sab^ ^bcy ^abcy ^ady ^bdy Scd)Pbc^p{x, Sbc)i (5.23) 

where, 

C ( X _ 2(s6d + Scd) / {Sab + ZSbc)\ 

^d;abc[Z, Sab, Sbc, Sabc, Sad, Sbd, Scd) — I ^ H I • IC).Z4J 

SabSad \ Sabc J 



A similar result holds for gluon c becoming soft, 

|^(...,a,6,c,...,e,...)|2-.P,,,;,|^(...,P,...,e,...)|^ (5.25) 

where, 

Pabc;e = Pd;abc{a ^ C,d^ e). (5.26) 

In the case that b is soft the matrix elements do not possess sufficient singularities since the 
collinear particles a and c are not directly colour-connected, 

|^(. . . , a, 6, c, . . .) 1^ less singular. (5-27) 

Here, there may be two powers of the small invariants in the denominator, but, when inte- 
grated over the appropriate (small) region of phase space this yields a vanishing contribution. 
On the other hand, for QED-like subamplitudes where the gluons are symmetrised there is a 
non-zero contribution when either gluon is soft. Note that in all cases where the soft particle 
is a quark or antiquark, there is also no singular contribution. 
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5.2.1 Soft/coUinecir limit of e'^e 5 partons 

For the specific case of the two-quark currents, the sum over all soft/coUinear limits is easily 
obtained, 



S,{Q,;1,2,3;Q,)V^^ 
S,iQ,;l,2,~3;Q,)V^ 



Pqi12;3 + ^Qi;123 + -Pl23;Q2 ^ "^l: 
(-Pqi;12Q2 + -^lQi3;Q2 ^SQil;2 
+-fQi;3Q22 + -^Qil2;Q2 -^3Q22;l) 
(A(3i2;Q2 + -^(3i;2Q23) 

P(l,2,3) 



23Q. 



(5.28) 



The only non- vanishing contributions in the soft / collinear limit from the four-quark currents 
are. 



2 / 

'I'lxiQu Q2; Q3, Qa] '^)y^ [^QilQiPQ^Qi^G + ^QilQiPqiQi^G 



(5.29) 



where the second term also includes both the unconnected soft /collinear contribution and 
interferences amongst the various subamplitudes. This is akin to the case of single soft gluon 
emission where we have. 



^(. ..,a,^,6, ...) 



^Pa-Pg Pb-Pg, 



A{...,a,b,...), 



and therefore. 



3fJ {A{. ..,a,g,b,...,c,d,...)) {A{. . . ,a,b, . . . ,c, g,d, . . .))^ 
Here, the soft factors are generated by the interference of the two eikonal factors. 

\Pa -Pg Pb-Pg) \Pc ' Pg Pd ' Pg J ' 



Adding up these limits for the five parton squared matrix elements gives. 



1 
3! 



+ 
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X 



-Pqi12;3 + -Pqi;123 + P] 



123;Q2 1;23Q 



+ -Pqi12;Q2 + -^3Qil;2 + P3Q^2;l) + a^4 (Aqi2;Q2 + -Pqi;2Q23) 



jY2 V Qi;i2Q2 ^ Qii2;Q2 ' " 3gii;:i ' -3g22;iy ' js^i 



ATS 



IQ4Q3 -^Qi;iQ3Q4 -^Q2;iQ4Q3 + -^Q2;1Q3Q4 



.(5.30) 



5.3 Two soft gluons 



Finally, we consider the contributions where two colour connected gluons are simultaneously 
soft. This was first studied by Berends and Giele [p5l and we include this contribution here 



for the sake of completeness. Similar results have been discussed by Catani [37]. For gluons 
h and c soft the colour ordered subamplitudes factorise, 



|^(. . . , a, 6, C, . . . , ) 1^ Sabcd{Sad, Sab, Scd, She, Sabc, Sbcd) \ A{. . . , a, d, . . .) |^. 

where the connected double soft gluon function is given by, 



(5.31) 



Sabcdi^ad, Sab, Scd, Sbc, Sabc, Sbcd) 



8s^ 



ad 



+ 



8s 



ad 



Sbc 



+ 



+ 



S abS bcdS abcS cd 

1 4 ^ 8(1-6) 



SabScd SabSbcd S^dSabc SabcSbcd 



"be 



Sabc SfiQd 



Here a and d are the hard partons surrounding the soft pair and may either be gluons or 
quark/ ant iquarks. In four dimensions, the double soft factor can be extracted from p5| by 
squaring and summing the helicity amplitudes for two adjacent soft gluons. Alternatively, 
it can be obtained by explicitly taking the double soft limit of squared matrix elements for 
processes involving more than two gluons. 



5.3.1 Double soft limit of e+e 5 partons 



As before, for the specific case of the two-quark currents, the connected double soft limit is 
easily obtained and summing over all contributions we find. 



S,{Qi;l,~2,~3;Q,)V>^ 



('S'qi123 + '5' 

s. 



123Q2 

2 



/112Q2 ^^ 
less singular 



(5.33) 
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Combining these limits yields the double soft singular factor for the full squared matrix 
elements, 



1 
3! 



+ 



(5.34) 



6 Summary 



In this paper we have examined the factorisation properties of squared tree level matrix 
elements when two particles are unresolved within the framework of QCD. At next-to-next- 
to-leading order, this knowledge is required for the analytic isolation of infrared singularities 
of n + 2 parton scattering processes and subsequent numerical combination with the one- loop 
n + 1 parton (single unresolved particle) and two-loop n parton contributions. 

The unresolved particles may be either soft gluons or groups of collinear particles or 
combinations of both. There are four double unresolved cases; two soft gluons, three simul- 
taneously collinear particles, two independent pairs of collinear particles and one soft gluon 
together with a collinear pair. In section ^ we reviewed the (trivial) factorisation that occurs 
when the unresolved particles are colour "unconnected". Such factorisation is well known 
and involves only the familar eikonal and Altarelli-Parisi splitting kernels used to describe 
single unresolved emission (see sect. |). 

When the unresolved particles are all colour "connected" , we find a similar factorisation. 
In section |^ we introduced new functions to describe the triple collinear and soft / collinear 
limits in addition to recalling the known double soft gluon limits of Berends and Giele ||2^ . 
These functions are universal and apply to general multiparton scattering amplitudes. As a 
check on our results, we find that the triple collinear splitting functions obey an expected 
N = 1 SUSY identity. In addition, in the strong ordered limit, where one particle is much 
more unresolved than the other, these factors become simple products of single unresolved 
factors, one associated with each unresolved particle. 

To illustrate the use of these double unresolved approximations, we have examined the 
singular limits of the tree level matrix elements for e'^e'^ —>■ 5 partons 0. In each case, 
we find that in the singular limit, the matrix elements can be approximated by a singular 
factor multiplying the tree level e~^e~ 3 parton matrix elements. These approximations 
will be of use in evaluating the 0(af) corrections to the three jet rate in electron positron 
annihilation. To achieve this however, much work still remains to be carried out. One 
important ingredient is to analytically integrate the approximations over the unresolved 



regions of phase space. A first step in this direction has been carried out in ref. where 
the hybrid subtraction method of has been used to evaluate the double unresolved 
singular contributions associated with a photon-gluon-quark cluster. There is in principle 
no reason why this approach should not be extended to the more general cases discussed here. 
A further ingredient is the evaluation of the two loop e~^e~ 3 parton matrix elements. 



22 



This is a formidable task in its own right and requires analytic expressions for the two loop 
box graph with massless internal and external legs for arbitrary dimension. So far, no such 
expression has been found. 
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